Resilience of <i-wave superconductivity to nearest-neighbor repulsion 
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Many theoretical approaches find d-wave superconductivity in the prototypical one-band Hubbard 
model for high-temperature superconductors. At strong-coupling (U > W, where U is the on-site 
repulsion and W = 8t the bandwidth) pairing is controlled by the exchange energy J = 4t 2 /U. One 
may then surmise, ignoring retardation effects, that near-neighbor Coulomb repulsion V will destroy 
superconductivity when it becomes larger than J, a condition that is easily satisfied in cuprates for 
example. Using Cellular Dynamical Mean-Field theory with an exact diagonalization solver for the 
extended Hubbard model, we show that pairing at strong coupling is preserved, even when V 2> J, 
as long as V < (7/2. While at weak coupling V always reduces the spin fluctuations and hence 
d-wave pairing, at strong coupling, in the underdoped regime, the increase of J = 4t 2 /{U - V) 
caused by V increases binding at low frequency while the pair-breaking effect of V is pushed to high 
frequency. These two effects compensate in the underdoped regime, in the presence of a pseudogap. 
While the pseudogap competes with superconductivity, the proximity to the Mott transition that 
leads to the pseudogap, and retardation effects, protect d-wave superconductivity from V. 



I. INTRODUCTION 

The existence of d-wave superconductivity in the one- 
band two-dimensional Hubbard model has been estab- 
lished through a variety of theoretical methods at both 
wealP^^U an( j strong coupling^H^I, in other words for 
one-site interaction U either much smaller or much larger 
than the bandwidth W. Generalizations of Dynamical 
Mean-Field Theory are particularly suited for the strong 
coupling limit, but they are also an excellent guide to 
the physics at weak to intermediate coupling! 33 " 42 !, These 
calculations suggest that pairing is maximized at inter- 
mediate coupling, where the on-site interaction U is of 
order the bandwidth W = 8t. Some non-pcrturbative 
calculations based on weak coupling ideas even agree at 
intermediate coupling 4 ^ with strong-coupling based ap- 
proaches. 

In all these approaches, spin fluctuation s with e ither 
an antiferromagnetic or a singlet characte r^ 8 * 44 * 45 ! have 
been argued to drive the pairing. These spin fluctuations 
result from the presence of an on-site Coulomb repulsion 
U. At strong coupling, the characteristic energy scale of 
these fluctuations, the exchange interaction J, is given 
by 4t 2 /U, and the d-wave gap symmetry adopted by the 
Cooper pairs allows them to avoid the direct effect of the 
on-site repulsion U . 

Little attention has been paid so far to the effect of the 
nearest-neighbor Coulomb repulsion (or extended Hub- 
bard interaction) V that we expect to be detrimental 
to d-wave superconductivity. Roughly speaking, in a 
simple BCS picture that does not take retardation into 
account, we expect the effective interaction to be the 
difference J — V. In ordinary phonon-mediated super- 
conductivity, the repulsion V is replaced by a smaller 
pseudopotential V c to account for the fact that bind- 
ing occurs at low frequencies through phonons while 
the Coulomb interaction acts over a broa d ener gy scale. 
This so-called Anderson-Morel mechanis m 146 * 47 !, leads to 



the following estimate for the Coulomb pseudopotential 
V c = V/(l + N(0)V\n(E F /uj D )) where the Debye fre- 
quency is ujd and the Fermi energy Ep. One expects 
that in strongly-correlated superconductivity, the ratio 
Ef/ojd must be replaced by a number closer to unity 
in which case this mechanism would no-longer be ef- 
fective and superconductivity should disappear as soon 
as V > J. This issue is crucial to understand high- 
temperature superconductors since that condition, or the 
weaker condition V > A s with A s the spin gap, is likely 
to be satisfied in these materials. From the value of the 
near-neighbor Coulomb interaction with a relative dielec- 
tric constant of order 10 we estimate V ps 400 meV while 
J fa f 30 me^. 

So far, it has been shown using a variational wave- 
function approach for the t — U — J — V model at strong 
coupling that superconductivity persists as long as^SI 
3 J/4 > S 2 (V- J/4) where 5 is the doping. A large-N cal- 
culation gives superconductivity at least up to V — 2J^ 
while Density Matrix Renormalization Group (DMRG) 
calculations on Hubbard or t — J laddera^"^ suggest 
that pairing can survive up to V ps 4 Jp^ At very small 
coupling (U <C W), it has been arguecP^that pairing is 
destroyed as soon as V > U(U/W). This weak coupling 
bound is close to the result of a FLEX calculation.^ 

Here we show that d-wave superconductivity in the 
one-band two-dimensional Hubbard model at strong cou- 
pling is in fact more robust than expected. Even for 
V ^> J, as long as the inequality V < U/2\s satisfied, su- 
perconductivity persists. This illustrates differences be- 
tween pairing at weak and at strong coupling, especially 
in the pseudogap regime. The resilience of d-wave super- 
conductivity to V can be traced to the increase in the 
effective J caused by V at strong coupling, i.e. when at 
half-filling the system is a Mott insulator. This increase 
in J is visible in the dynamics of both the spin suscep- 
tibility and the Cooper pair Green's function which is 
enhanced by V at low frequency. The pair-breaking ef- 
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feet of V manifests itself at higher frequency. This leads 
overall to a sizeable range of values of V where the order 
parameter is essentially independent of V in the under- 
doped regime where a pseudogap appears. Recall that 
at strong coupling, the pseudogap extends up to optimal 
doping, a sign that Mott physics extends well beyond 
half-filling^. 

We use Cellular-Dynamical Mean-Field Theory that 
allows us to study the strong-coupling limit and allows 
one to take into account both J and V as well as the 
effect of retardation. 

We first present the method and model and conclude 
with a discussion after the presentation of the results 
obtained from large scale numerical calculations. 
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FIG. 1. Cluster and bath parametrization used in this work. 
Hopping terms are shown by full lines, bath hybridization by 
dashed lines and bath pairing terms by dotted lines. 



II. MODEL AND METHOD 

We start from the one-band extended Hubbard model 
on a square lattice, 

where cfl is the destruction (creation) operator for an 
electron of spin a at site i and n ia — cj a c i a . is the 
corresponding number operator (n, = + n^). We 
assume a band structure close to that of YBa2Cu307, 
with nearest-neighbor hopping t set to unity, diagonal 
hopping t' — —0.3 and third-neighbor hopping t" = 0.2, 
unless otherwise indicated. This model is solved with 
a 4 site plaquette Cellular Dynamical Mean Field The- 
ory (CDMFT) at T — using an exact diagonalization 
solveJ32E2HHl that allows us to obtain real-time quantities 
without analytic continuation. This approach has been 
used to reveal the presence of d-wave superconductivity 
in the one-band Hubbard modelp^l and to study the pair- 
ing dynamics and retardation when V = The results 
obtained for d-wave superconductivity with the plaque- 
tte are essentially identical to those obtained with larger 
clusters.^ 

As in previous studies of the 2D Hubbard model with 
this approach^, the plaquette is hybridized with a set of 
8 bath orbitals, as illustrated in Fig.[T] with 6 parameters 
to be determined through the CDMFT self-consistency 
relation: the hybridizations #i : 2 between the bath and the 
cluster, the bath energies £i : 2 and the bath d-wave pair- 
ing parameters which, when nonzero, signal the pres- 
ence of superconductivity and lead to a non-vanishing 
cluster d-wave order parameter 

tp = c ^ c 3i - c A°ji + c - c - ( 2 ) 

The order parameter is extracted from the lattice 
Green's function while the doping is measured on the 
cluster. The distance function is defined in Ref. [311 
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FIG. 2. (Color online) Blue dots: Inverse frequency of the 
dominant peak in the antiferromagnetic susceptibility of the 
extended Hubbard model at half-filling, as a function of the 
nearest-neighbor interaction V, for U = 16 and U — 32 (t = 
1). These inverse frequencies are expected to be J" 1 = (U — 
V)/4t 2 in the large-t/ limit. 



We use cutoff u c — 2 and fictitious inverse temperature 
\3 = 50. 

While the extended interaction V within the cluster 
can be treated exactly with an approach like CDMFT, 
the coupling to neighboring clusters requires the Hartree 
approximation.^! More specifically, the model Hamilto- 
nian Eq. ([!]) is modified by replacing the inter-cluster 
interaction by 

V riinj + V {fiiTij + mfij — fiifij) (3) 

{i,j)c <M>ic 

where (i,j) c denotes nearest-neighbor pairs within the 
plaquette and (i,j)i c nearest-neighbor pairs across pla- 
quettes. The mean field fii must be determined self- 
consistently; in practice, it is treated like the six bath 
parameters wihtin the CDMFT self-consistency loop and 
thus both the dynamical mean field (represented by the 
bath parameters) and the static Hartree mean field are 
converged simultaneously. 

To verify the accuracy of this approach, consider 
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FIG. 3. (Color online) d-wave order parameter ip obtained 
from the off-diagonal component of the lattice Green's func- 
tion as a function of cluster doping for U — 4, 8 and 16 and 
various values of V. 



momentarily the half-filled, particle-hole symmetric ex- 
tended Hubbard model with t' = t" = 0. In that simple 
case, the Hartree field fii is fixed by particle-hole sym- 
metry. We find that the dynamical spin susceptibility at 
the antiferromagnetic wavevector Q = (n, tt) has a dom- 
inant peak at a low frequency uj x . The position of this 
peak is shown on Fig. [5] as a function of V for two values 
of U. In the large U limit, the half-filled extended Hub- 
bard model should map to the Heisenberg model with a 
super-exchange parameter J = At 2 /(U — V). The energy 
denominator is easily understood by comparing the en- 
ergy, at t = 0, between a configuration where all sites 
are exactly singly-occupied and another configuration in 
which one electron has vacated a site in order to doubly 
occupy a neighboring site. Four V bonds are lost, but 
three are gained. Since we expect the frequency lu x to 
scale like J in the large-?/ limit, Fig. [2] confirms that we 
obtain the correct scaling with V. Hence, our approach 
leads to the correct strong-coupling physics. 



III. RESULTS 

To concentrate on the effect of short-range spin fluctu- 
ations on superconductivity, we ignore the possibility of 
long-range antiferromagnetic order predicted by CDMFT 



at low dopingpS as well as the possibility of long-range 
charge density waves in the presence of V . At weak cou- 
pling charge-de nsity w aves may occur for U ~ 4V in 
two-dimensionsP 5 * 61 * 62 ! We found a similar result at half- 
filling for strong coupling. 

The (i-wave order parameter for YBCO hopping pa- 
rameters is displayed on Fig. [3] for hole-doping at U = 4, 
8 and 16. In each panel, obtained for a given value of 
U, we display the results for different values of V. The 
Mott transition at half-filling, that separates weak from 
strong coupling, occurs around U = QpSiMi As observed 
before,^! in the weak coupling case (U = 4) superconduc- 
tivity is strongest at half-filling when competition with 
antiferromagnetism is prohibited (Fig. [3^,). At larger cou- 
pling, the Mott transition destroys superconductivity at 
half-filling (Figs.|]b,c). 

The value of V influences the order parameter in strik- 
ingly different ways at weak and at strong coupling. For 
[7 = 4, superconductivity has essentially disappeared at 
V/U = 1.5/4 = 0.375, in agreement with the upper limit 
V/U = U/W found by weak coupling analysis.^ This 
is consistent with the fact that at weak coupli ng, V al- 
ways decreases the strength of spin fluctuations! 61 * 62 ! By 
contrast, at strong coupling, V can increase the strength 
of spin fluctuations through J = 4t 2 /(U — V), and for 
the same ratio V/U — 0.375, the order parameter is still 
large for U = 8 on Fig. [3]d . In fact, it is barely influ- 
enced by V in the pseudogap (underdoped) regime close 
to half-filling (x — 0.05). On the lower panel, for U = 16, 
we have V/ J = 16 at V = 8 and superconductivity still 
persists. The ratio V/J is maximum at V = U/2. We 
estimate that V « U/2 is the upper bound for the the 
appearance of <i-wave superconductivity (ignoring com- 
peting orders). 

The link between superconducting order parameter 
and spin fluctuations was demonstrated at V = by the 
correlation between the peaks in the imaginary part of 
the an omalo us self-energy and the peaks in the spin fluc- 
tuatiorPSES This correlation persists here (not shown). 
To further investigate this link between spin fluctua- 
tions and superconductivity, let us thus first focus on 
the spin dynamics revealed by the the local spin spectral 
function x"(w) illustrated on Fig. [4] Panel (A) shows 
x"(cj) at U = 8, V = 1 for the underdoped (x = 0.05) 
and overdoped (x = 0.20) regimes. The corresponding 
charge susceptibility is negligible on this scale. Spin 
fluctuations are much smaller in the overdoped regime 
and their weight is spread in wavevector contrary to the 
underdoped case. Panel (B) shows that in the under- 
doped regime [x = 0.05) where there is a pseudogap, 
a moderate extended interaction V at strong coupling 
mainly shifts the spectrum of Q = (n, tt) spin fluctua- 
tions to higher frequencies, as expected from the increase 
of J = At 2 / (U — V) caused by V; it decreases somewhat 
the overall spectral weight without affecting the spectral 
shape. 

To show that low-frequency spin fluctuations are ap- 
parent in the pair dynamics, we compare the character- 
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FIG. 4. (Color online) (A) Imaginary part of the local spin 
susceptibility x"( w ) f° r U = 8, V 7 = 1 for an underdoped 
(x — 0.05) and an overdoped (x = 0.20) case. (B) In the 
underdoped regime (x = 0.05) apart from a shift of the low- 
frequency peak to higher frequencies, it is mostly the ampli- 
tude, not the structure of the imaginary part of the plaquette 
antiferromagnetic spin susceptibility x"(Q,^) that is slightly 
affected by V . 



istic frequency w x of spin fluctuations, given by the po- 
sition of the dominant peak of x"(Q,w), with a charac- 
teristic frequency in the pair dynamics, as was done in 
Ref. 1301 The pair dynamics can be studied through the 
integral 



du>' 

7T 



(4) 



where F R is the retarded Gork'ov function (or anoma- 
lous Green's function) defined in imaginary time by 
Fij = — (Tcif(r)cj±(0)) with i and j nearest-neighbors. 
The infinite frequency limit of If{oj) is equal to (c^Cj^) 
which in turn is proportional to the T = ci-wave order 
parameter tp. If(u>) is useful to estimate the frequencies 
relevant for binding. We call If{w) the cumulative order 
parameter. 

In BCS theory, If(uj) is a monotonically increasing 
function of cu that reaches its asymptotic value at the 
BCS cutoff frequency uj c W^ In the Eliashberg approach 
that includes retardation as well as the Coulomb pseu- 
dopotentialp^ the function overshoots its asymptotic 
value at frequencies near the main phonon frequencies 
before decaying to its final value because of pair breaking 
effects at higher frequencies. We define the characteris- 
tic frequency top as the point where If(w) reaches half 
of its asymptotic value. If one imagines that F^(u) is 
made of a single peak, then that peak would be located 
at up . Figure [5] illustrates the cumulative order param- 
eter If(uj) for various values of V. The top two panels 
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FIG. 5. (Color online) Integral of the anomalous Green's func- 
tion (or Gork'ov function) Ip(cj) obtained after extrapolation 
to r\ — of uj + ir\ for several values of V at (A) U = 8, 
x = 0.05 (B) U = 8, x = 0.2 (C) U = 16, x = 0.05. The 
asymptotic value of the integral, 7f(oo), equal to the order 
parameter, is shown as horizontal lines. We call If{oj) the 
cumulative order parameter. The characteristic frequency ujf 
is defined as the frequency at which If(oj) is equal to half 
of its asymptotic value. The horizontal arrow in panel (A) 
indicates how ujf is obtained. 



are for two values of doping x at U — 8, and the bottom 
one for U — 16 in the underdoped regime. The asymp- 
totic value is indicated by a horizontal line. Following a 
sharp rise around uj — up, the function has a maximum 
and then decreases towards its asymptotic value, as in 
Eliashberg theory. 

The link between the characteristic spin-fluctuation 
frequency lo x and the characteristic frequency of the pair 
up is summarized on Fig. |6j which shows the evolu- 
tion with doping of: (A) the position uj x of the main 
AF susceptibility peak, (B) its strength, and finally (C) 
the characteristic frequency lof, for four values of V at 
U = 8. Clearly, uj x , \ an d wj? all decrease with doping. 
This shows that spin fluctuations and pair dynamics are 
strongly linked. 

Even for V — however, the increase in J as mea- 
sured by u> x and up in the underdoped regime, x < 0.10, 
does not correspond to an increase in the order parame- 
ter. As can be seen in Fig. [3j for x = 0.05 and x = 0.20 
the d-wsve order parameter is essentially identical when 
V = and the increase in ujp near half- filling manifests 
itself only in the maximum value of the cumulative or- 
der parameter Ip(oj) in Fig. [5j not in the value of the 
order parameter. Clearly the order parameter does not 
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FIG. 6. (Color online) Panel (A) position of the lowest peak 
in the AF susceptibility x"^) an d t ne strength \ of that peak 
on Panel (B), for (7 = 8. Panel (C) shows the corresponding 
SC characteristic frequency uif- 



increase in the underdoped regime despite the increase in 
the characteristic spin fluctuation frequency lu x because 
of another effect. That effect is the pseudogap. The pseu- 
dogap removes many of the states near the Fermi energy 
that would otherwise be paired. 

The presence of the pseudogap can be seen by con- 
trasting the single-particle local density of states at two 
dopings. In the overdoped regime (x = 0.20) the local 
density of states has no pseudogap in the normal state, 
as seen from the dashed red line on Fig.[7)3. Only the su- 
perconducting gap is visible. By contrast, even at V = 
there is clearly a pseudogap in the underdoped regime, as 
can be seen from Fig. [7^,. The effect of superconductivity 
manifests itself only through the small more symmetric 
gap near uj = 0.^ The pseudogap, as measured from the 
peak to peak distance, increases slightly with V in the 
underdoped regime. 

The detrimental effect of the pseudogap on the super- 
conducting order parameter was demonstrated in Fig. 7 
of Ref. 1301 where the bubble contribution to the pair- 
ing susceptibility decreases as we approach half-filling. 
This means that in the absence of interactions between 
the particles forming the pair (represented by vertex cor- 
rections), self-energy effects disfavor superconductivity. 
This physics is also present at weak to intermediate cou- 
pling when the pseudogap is induced by long wavelength 
antiferromagnetic fluctuations^ (and not by the incipi- 
ent shor t-ran ge Mott localisation that appears at strong 
couplings!). 



IV. DISCUSSION 
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FIG. 7. (Color online) Local density of states for four values 
V = 1,2, 3,4 and U = 8 for (A) underdoped x = 0.05, and 
(B) overdoped x = 0.2 regimes. The dotted line in (B) is for 
the normal state V = 0, all other lines in (A) and (B) are for 
the superconducting state. 



The effect of spin fluctuations and of near-neighbor re- 
pulsion V on pairing is relatively straightforward in the 
overdoped regime (a; > 0.10). There, as seen on Fig. ^ijp 
and Figs. |6^i,c, the order parameter increases with the 
characteristic spin frequency lo x and the characteristic 
pair frequency cup. Like for the V — case, superconduc- 
tivity disappears at finite V when the amplitude x °f the 
lowest frequency peak in the spin fluctuations (Fig. [6]d) 
vanishes,^ as observed in experiment 00323 The main ef- 
fect of V is to decrease the spin fluctu ations, as happens 
systematically at weak coupling! 61 * 62 ! 

As we approach half-filling, the decrease of the order 
parameter for all values of V (Figs. [3j},c), despite the 
increase in the amplitude x °f the lowest peak in the 
spin fluctuations (Fig. [6]), is mainly due to the increase 
of the pseudogap: it removes more and more states from 
the Fermi energy as the doping x decreases (Fig. [7| . Note 
that V causes opposite changes in the amplitude x an d 
in the characteristic frequencies (Figs [6]), which might 
explain the near independence of the order parameter 
with respect to V in Figs. [3j3,c. However, these changes 
are small on a relative scale and are probably not the 
main reason for the insensitivity of the order parameter 
to V. 

The resilience of d-wave superconductivity to V at 
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strong coupling is best understood from the U = 16 re- 
sults for the cumulative order parameter Ip (oj) shown in 
Fig. [5b for x = 0.05. The largest value of If(u>) scales 
roughly like J — At 2 / (U — V), in other words V increases 
the binding at low frequency, where the retardation is 
large. However, V also has the expected pair-breaking 
effect: larger V causes larger pair-breaking so that the 
low frequency increase in If(uj) is essentially compen- 
sated by the time If(u>) reaches its asymptotic value (the 
order parameter). 

To understand more deeply the dual role of V, as both 
pair binding and pair breaking, it is helpful to return to 
the solution of the Eliashberg equations in the electron- 
phonon case. Figs. 7, and 13 of Ref. |47] show that the 
imaginary part of the gap function remains positive for 
all frequencies when the Coulomb pseudopotential van- 
ishes. By contrast, one verifies from Fig. 11 of the same 
paper that the imaginary part of the gap function changes 
from positive to negative as frequency increases when the 
Coulomb pseudopotential is finite. A sign change in the 
imaginary part of the gap function should lead to a sim- 
ilar behavior in the imaginary part of the Gorkov func- 
tion F(uj). In this case its integral -Tf(w), as defined in 
Eq. Q, should reach a maximum before decreasing to 
its asymptotic large frequency value. By contrast, when 
the imaginary part of the gap function remains positive 
at all frequencies, i.e. in the absence of Coulomb pseu- 
dopotential, If{uj) should increase monotonically to its 
asymptotic value, as in BCS theory. Since we observe 
that If(ui) has a maximum at finite frequency, we are 
in the case where the Coulomb pseudopotential is impor- 
tant. Comparison of Figs. 11 and 13 of Ref. ITfl also shows 
that the maximum value of the imaginary part of the gap 
function (and that of its frequency integral, as one can 
check) increases with the electron-phonon coupling con- 
stant. In our case, the increase in J = 4t 2 /(U — V) 
caused by V similarly leads to an increase in the maxi- 
mum value of the cumulative order parameter Ip{u). In 
other words, we can surmise that the increase in V leads 
to larger binding at low frequency, like an increase in the 



electron-phonon coupling constant, but it also leads to 
more pair-breaking at large frequency, coming from the 
Coulomb pseudopotential-likc effect of V, and these two 
effects nearly cancel each other. 

In the electron-phonon case, one also notes that the 
imaginary part of the gap function vanishes at about two 
to three times the maximum phonon frequency. This 
is where the cumulative order parameter If{uj) would 
reach its asymptotic value. In the case of spin fluctua- 
tions, 2 J is a measure of the width of the spin-fluctuation 
spectrum. Comparing Fig. [5^, for U = 8 and Fig. [5b 
for U — 16, one notices that the frequency for which 
If(lu) reaches its asymptotic value is smaller by a factor 
of roughly two for U = 16 where the V = value of J is 
twice as small as for U = 8. 

The opposing effects discussed above conspire to leave 
the order parameter rather insensitive to V in the pres- 
ence of a pseudogap. The pseudogap induced by Mott 
physics^ is detrimental to superconductivity, but in its 
presence superconductivity is effectively protected from 
near-neighbour repulsion V. Indeed, in a doped Mott 
insulator, short-range incipient localisation is strong 
enough to create a pseudogap, but while V causes pair- 
breaking at high frequency, it also enhances spin fluctu- 
ations (J = At 2 /(U — V)) at low frequencies, thus com- 
pensating the pair- breaking effect. Overall, retardation 
effects are crucial for the resilience of ci-wave supercon- 
ductivity to near-neighbor repulsion. 
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